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/"*\ , We propose a new backbone-torsion-energy term in the force field for protein systems. This 

torsion-energy term is represented by a double Fourier series in two variables, the backbone di- 
CNj ' hedral angles (f> and ij). It gives a natural representation of the torsion energy in the Ramachan- 

04 , dran space in the sense that any two-dimensional energy surface periodic in both cf> and ip can 

be expanded by the double Fourier series. We can then easily control secondary-structure-forming 
tendencies by modifying the torsion-energy surface. For instance, we can increase/decrease the 
a- helix- forming-tendencies by lowering/raising the torsion-energy surface in the a- helix region and 
. likewise increase/decrease the /3-sheet- forming tendencies by lowering/raising the surface in the /3- 

sheet region in the Ramachandran space. We applied our approach to AMBER parm94 and AMBER 
parm96 force fields and demonstrated that our modifications of the torsion-energy terms resulted in 
the expected changes of secondary-structure-forming-tendencies by performing folding simulations 
of a-helical and /3-hairpin peptides. 
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I. INTRODUCTION 



e 

A force field (or potential energy) for protein systems is necessary to perform molecular simulations based on 
Monte Carlo (MC) and molecular dynamics (MD) methods. Well-known force fields are, for instance, AMBER 
II CHARMM ||, OPLS 00, GROMOS 0, and ECEPP 0. These force fields have been parameterized to 
fit experimental data of small molecules and, for some terms, the results of quantum chemistry calculations and MD 
simulations. We have recently carried out detailed comparisons of three versions of AMBER (parm94 0, parm96 
0, and parm99 0), CHARMM 0, OPLS-AA/L 0, and GROMOS by generalized-ensemble simulations of 
two small peptides in explicit solvent. 0, The results of these comparisons indicated that these force fields 
yield quite different secondary-structure-forming tendencies. It was shown that AMBER parm94 is the most (and 
too much) a-helix-forming among the six force fields studied and that AMBER parm99 and CHARMM give ample 
amount of a-helix structures, whereas AMBER parm96, OPLS-AA/L, and GROMOS are more /3-sheet-forming than 
the rest. l lCtllll These results were confirmed by the folding simulations of the two peptides with implicit solvent 
I/"") ' model. |l2Ul3lll4| Our conclusion was that we need to refine and improve the existing force-field parameters in order 
to yield the secondary-structure-forming tendencies that agree with experimental implications. jToL ITlL IT3. lilj. Il4| 
Among the force-field terms, the torsion-energy term is the most problematic. For instance, the parm94, parm96, 
and parm99 versions of AMBER differ mainly in the backbone-torsion-energy parameters. The main changes from 
OPLS-AA to OPLS-AA/L can also be found in the torsion-energy term. Other trials of force-field refinement 
mainly concentrate on the torsion-e nergy t erms. These modifications of the torsion energy are usually based on 
quantum chemistry calculations |l5l ll& Il7l Hsl Il9| . It was also proposed to set the backbone torsion-energy term 
simply to zero. [2(j We have proposed another force-field refinement procedures 0,0,0]- Using the standard 
functional forms of the force fields, we refine the force- field parameters so that they are the most consistent with 
the coordinates of proteins in the Protein Data Bank (PDB) database. This is achieved by minimization of forces 
• i-H , acting on the atoms with the coordinates of the PDB with respect to the force-field parameters. We found that our 
refinement of partial-charge and backbone-torsion-energy parameters resulted in the improvement of three versions of 
AMBER (parm94, parm96, and parm99), CHARMM, and OPLS-AA [H[llQ3. 

In this article, we propose a new backbone-torsion-energy term, which is represented by a double Fourier series in 
two variables, the backbone dihedral angles <p and ip. This expression gives a natural representation of the torsion 
energy in the Ramachandran space pH in the sense that any two-dimensional energy surface periodic in both <fi and tjj 
can be expanded by the double Fourier series. We can then easily control secondary-structure-forming tendencies by 
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modifying the backbone-torsion-energy surface. We accommodated our approach to AMBER parm94 and AMBER 
parm96 and tested whether our modifications of the backbone-torsion-energy term resulted in the expected change of 
secondary-structure-forming-tendencies by performing folding simulations of a-helical and /3-hairpin peptides. 

In section 2 the details of the double Fourier series as the backbone-torsion-energy term are given. In section 3 the 
results of applications of the double Fourier series to modifications of AMBER parm94 and parm96 force fields are 
presented. Section 4 is devoted to conclusions. 



II. METHODS 



The existing all-atom force fields for protein systems such as AMBER and CHARMM use essentially the same 
functional forms for the conformational potential energy E con f except for minor differences. E con { can be written as, 
for instance, 

Econf = E-BL + EbA + -Etorsion + E/nonbond ■ (1) 

Here, -EbLj -E^ba, £- torS i 0n , an d -E^onbond represent the bond-stretching term, the bond-bending term, the torsion-energy 
term, and the nonbonded energy term, respectively. The torsion energy is usually given by 

^torsion = Ef[ 1 + C0S ( n$ ~ 7n)] , (2) 

dihedral angle <E> n 

where the first summation is taken over all dihedral angles <J> (both in the backbone and in the side chains), n is 
the number of waves, 7„ is the phase, and V n is the Fourier coefficient. Separating the contributions E(cp,ip) of the 
backbone dihedral angles and ip from the rest of the torsion terms E rest , we can write 

EVsion = E(4>, ip) + E rcst , (3) 

where we have 

S(0,^)=^^[l + co S (m0- 7 m)]+^Y [1 + COs(n ^~ 7n)] ■ (4) 

m n 

For example, the coefficients for the cases of AMBER parm94 and AMBER parm96 are summarized in Table [I] and 
we can rewrite E(<p : ip) in Eq. (0} as follows: 

E paim94 ,(<p,ip) = 2.7- 0.2 cos 20 -0.75 cos tp - 1.35 cos 2ip - 0.4 cos 4^ , (5) 



£Wm96(0,-0) = 2.3 + 0.85 cos 0-0.3 cos 20 + 0.85 cos tp-0.3 cos 2ip . (6) 

The backbone-torsion-energy term E(4>, tp) in Eq. Q is a sum of two one-dimensional Fourier series: one is for and 
the other for tp. The two variables and ip are independent, and no correlation of and ip can be incorporated. Any 
periodic function of and ip with period 2tt can be expanded by a double Fourier series. As a simple generalization 
of E(<j),ip), we therefore propose to express this backbone torsion energy by the following double Fourier series: 



£(<fi,ip) = a + (b m cos m0 + c m sin mq 

m— 1 

oo 

+ (d n cos mp + e n sin nip) 



71=1 

CO CO 



+ (fmn cos m<p cos 7i ip + g mn cos m<p sin nip 

m—l n—1 

+ h mn sin m0 cos nip + i mn sin rrup sin nip) . (7) 

Here, m and n are the numbers of waves, a, b m , c m , d n , e n , f mn , 9mn, h mn , and i mn are the Fourier coefficients. This 
equation includes cross terms in and ip, while the original term in Eq. Q has no mixing of and ip. Therefore, 
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our new torsion-energy term can represent more complex energy surface than the conventional one. The Fourier 
coefficients, by definition, are given by 

-i />7T />7T 2 1 /"180 plSO 

where a are the normalization constants and x(cj>, ip) are the basis functions for the Fourier series. TablelTllsummarizes 
these coefficients and functions. Here, <p and ip are given in radians, and <j6 and ip are in degrees (<p = jgQ<f>, ip — ifgVO- 
Hereafter, angular quantities without tilde and with tilde are in radians and in degrees, respectively. 

Finally, £(<p, ip) in Eq. Q and -E rcs t in Eq. © define our torsion-energy term in Eq. (instead of Eq. J2J): 

^torsion = £{4>, Ip) + E Iest . (9) 

The double Fourier series in Eq. (JJJ) is particularly useful, because it describes the backbone-torsion-energy surface 
in the Ramachandran space. The Fourier series can express the torsion-energy surface £(<p,ip) that was obtained by 
any method including quantum chemistry calculations. 0, 0, 0, 0, [l8L IT^ | 

Moreover, one can refine the existing backbone-torsion-energy term and control the secondary-structure-forming 
tendencies of the force fields. For example, a-helix is obtained for ((p,ip) w (—57°,— 47°), 3io-helix for (4>,ip) ~ 
(-49°, -26°), 7r-helix for (0,^) » (-57°, -70°), parallel /3-sheet for (</>>) ~ (-119°, 113°), antiparallel /3-sheet for 
((pj'ip) ~ (—139°, 135°), and so on. [2l| Hence, if the existing force field gives, say, too little a-helix-forming tendency 
compared to experimental results, one can lower the backbone-torsion-energy surface near ((p,if>) = (—57°,— 47°) in 
order to enhance a-helix formations. 

We can thus write 



£{<P,i>) = E{<P,iP) - f{<P,iP) , (10) 

where E(cj>, ip) is the existing backbone-torsion-energy term that we want to refine and f(<j>, ip) is a function that has 
peaks around the corresponding regions where specific secondary structures are to be enhanced. There are many 
possible choices for f((p,ip). For instance, one can use the following function when one wants to lower the torsion- 
energy surface in a single region near (</>, ip) = (<p , ip ): 

[ , otherwise , 

where A, B, and ro are constants that we adjust for refinement. In this case, the energy surface is lowered by /(</>, ip) 
in a circular region of radius r , which is centered at (<p,ip) — (4>o,ipo)- Note that we should also impose periodic 
boundary conditions on f{(p,ip). 

We then express £(<p,ip) in Eq. (|10|l in terms of the double Fourier series in Eq. @, where the Fourier coefficients 
are obtained from Eq. (jHJ. Hence, we can fine-tune the backbone-torsion-energy term by the above procedure so that 
it yields correct secondary-structure-forming tendencies. 



III. RESULTS AND DISCUSSION 

We now present the results of the applications of our backbone torsion energy in Eq. (JJJ. In this section, we consider 
the following truncated Fourier series: 

£(4>,ip)—a + b\ cos 4> + c\ sin <f> + b% cos 20 + ci sin 1<p 

+ g?i cos ip + e\ sin ip + c? 2 cos 2ip + e 2 sin lip 

+ /n cos 4> cos ip + gu cos ^ sin ^ 

+ hu sin (p cos ip + in sin <f> sin ip . (12) 

This function has 13 Fourier-coefficient parameters. We will see below that this number of Fourier terms is sufficient 
for most of our purposes, but that for some cases more number of Fourier terms are preferred. 

We first check how well the truncated Fourier series in Eq. (|12|) can reproduce the original AMBER parm94 and 
AMBER parm96 backbone-torsion-energy terms in Eqs. JSJ and iJBJ. Because these functions are already the sum of 
one-dimensional Fourier series and subsets of the double Fourier series in Eq. (J7J , the Fourier coefficients in Eq. (JSJ can 
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be analytically calculated and agree with those in Eqs. (0 and JBJ except for the last one (that for cos 40) in Eq. J5J. 
This term is missing in Eq. (|12[) . These cases thus give us good test of numerical integrations in Eq. (JSJ. The numerical 
integrations were evaluated as follows. We divided the Ramachandran space (—180° < <j) < 180°, —180° < i/j < 180°) 
into unit square cells of side length e (in degrees). Hence, there are (360/e) 2 unit cells altogether. The double integral 

on the right-hand side of Eq. (JSJ was approximated by the sum of £ (^jgo4>, ifp^) x ( ifo^' ifo 1 ^) x (^) 2 ' wnere each 

^ (lfo^' ifo^) x ( llo^' lfo^J was eva l ua ted at one of the four corners of each unit cell. We tried two values of e (1° 
and 10°). Both cases gave exact agreement of Fourier coefficients with the results of the analytical integrations for at 
least six digits (see Tables ITTTI and llVl below) . 

In Fig. we compare the backbone-torsion-energy surfaces of the original AMBER parm94 and AMBER parm96 
with those of the corresponding double Fourier series in Eq. 1)12(1 . Hereafter, the primed labels for figures such as 
(a') indicate that the results are those of the double Fourier series. As can be seen from Figs, db) and^b'), the 
backbone-torsion-energy surfaces are in complete agreement for AMBER parm96, whereas we see a little difference 
for AMBER parm94 between Figs. ^Ja) and^a'). As discussed above, this slight difference for AMBER parm94 
reflects the fact that the cos 40 term in Eq. (0 is missing in the truncated double Fourier series in Eq. i|12|) ■ 

We now consider the double Fourier series of non-trigonometric functions. The functions are those in Eqs. i|10|) 
and Qllfl. We try to fine-tune the original AMBER parm94 and AMBER parm96 force fields by subtracting f(<fi, tji) 
in Eq. i|ll|) from the original functions. The criterion for fine-tuning is, for instance, whether the refined force fields 
yield better agreement of the secondary-structure-forming tendencies with experimental implications than the original 
ones. For this we need good experimental data. Because the purpose of the present article is to test whether or not we 
can control the secondary-structure-forming tendencies, we simply consider extreme cases where we try to modify the 
existing force fields so that desired secondary structures may be obtained regardless of the tendencies of the original 
force fields. Note that the original AMBER parm94 and AMBER parm96 favor a-helix and /3-sheet, respectively. 

The function /(0, 0) in Eq. (|llf> reduces the value of E((j>, 0) in a circle of radius r$ with the center located at 
((f>o,t(jo). We used fo = 100° and B = 5,000 (degrees) 2 . The coefficient A is calculated by Eq. (|llfl from the other 
parameters f{4>o,' t Po)i ^o, and B. Namely, we have 

A = /(0 O) ^o)exp (^\ . (13) 

We used (0oi0o) — (—57°, —47°) and (f/'cbV'o) = (—130°, 125°) in order to enhance a-helix-forming tendency and 
/3-sheet-forming tendency, respectively. The central values f(4>o,ipo) that we used were 3.0 kcal/mol and 6.0 kcal/mol 
for enhancing a-helix and /3-sheet, respectively, in the case of AMBER parm94. They were both 3.0 kcal/mol in the 
case of AMBER parm96. We remark that the large value of /(0Oj0o), 6.0 kcal/mol, that was necessary to enhance 
/3-sheet in the case of AMBER parm94 reflects the fact that the original force field favors a-helix so much. 

In Fig.|21we show the backbone-torsion-energy surfaces modified according to Eq. (|10|) . We see that Eq. I|10|) reduces 
the torsion energy in the circular regions that correspond to the a-helix region (al and bl) and the /3-sheet region 
(a2 and b2). Hence, there are four cases: a-helix is enhanced from the original AMBER parm94 (al) and AMBER 
parm96 (bl), and /3-sheet is enhanced from the original AMBER parm94 (a2) and AMBER parm96 (b2). 

These modified backbone-torsion-energy functions were expanded by the truncated double Fourier series in Eq. i|12|) 
by evaluating the corresponding Fourier coefficients from Eq. JSJ. For the numerical integration we again tried two 
values of the bin size e (1° and 10°). The obtained Fourier coefficients are summarized in Table ITTTI in the case 
of AMBER parm94 and Table ITVl in the case of AMBER parm96. For comparisons, the Fourier coefficients of the 
original AMBER force fields (before modifications) are also listed. We see that the two choices of the bin size e gave 
essentially the same results (up to at least 3 digits). 

In Fig. |3 we show the backbone-torsion-energy surfaces represented by the truncated double Fourier series. Com- 
paring these with the original ones in Fig. [21 we find that the overall features of the energy surfaces are well reproduced 
by the Fourier series. If more accuracy is desired, we can simply increase the number of Fourier terms in the ex- 
pansion. As we see below, the present accuracy of the Fourier series was sufficient for the purpose of controlling the 
secondary-structure-forming tendencies towards a-helix or /3-sheet. 

We examined the effects of the above modifications of the backbone-torsion-energy terms in AMBER parm94 
and AMBER parm96 (towards specific secondary structures) by performing the folding simulations of two peptides, 
C-peptidc of ribonuclease A and the C-terminal fragment of the Bl domain of streptococcal protein G, which is 
sometimes referred to as G-peptide j^. The C-peptide has 13 residues and its amino-acid sequence is Lys-Glu-Thr- 
Ala-Ala-Ala-Lys-Phe-Glu-Arg-Gln-His-Met. This peptide has been extensively studied by experiments and is known 
to form an a-helix structure Because the charges at peptide termini are known to affect helix stability p3ll24]|. 
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we blocked the termini by a neutral COCH3- group and a neutral -NH2 group. The G-peptide has 16 residues and 
its amino-acid sequence is Gly-Glu-Trp-Thr-Tyr-Asp-Asp-Ala-Thr-Lys-Thr-Phe-Thr-Val-Thr-Glu. The termini were 
kept as the usual zwitter ionic states, following the experimental conditions p3 . l2RL I2r3 | . This peptide is known to 
form a /3-hairpin structure by experiments |22ll25l I26]] . 

Simulated annealing [27| MD simulations were performed for both peptides from fully extended initial conformations, 
where the four versions of the truncated double Fourier series (which were described in Tables IIIII and IIVI and in 
Fig. |2| were used for the backbone-torsion-energy terms of AMBER parm94 and AMBER parm96 force fields. For 
comparisons, the simulations with the original AMBER parm94 and parm96 force fields were also performed. The 
unit time step was set to 1.0 fs. Each simulation was carried out for 1 ns (hence, it consisted of 1,000,000 MD steps). 
The temperature during MD simulations was controlled by Berendsen's method 28] . For each run the temperature 
was decreased exponentially from 2,000 K to 250 K. As for solvent effects, we used the GB/SA model [2!J,|3(j. We 
modified and used the program package TINKER version 4.1 31] for all the simulations. For both peptides, these 
folding simulations were repeated 60 times with different sets of randomly generated initial velocities. 

In Table the numbers of obtained conformations (final conformations) with a-helix structures and /3-hairpin 
structures are listed for the case of C-peptide, which is known to form a-helix structures by experiments. We used 
DSSP 32] for the criterions of secondary-structure formations. We see that all 60 conformations obtained from the 
simulations using the original AMBER parm94 are a-helix structures and that for the original AMBER parm96, 14 
out of 60 conformations are a-helix structures and 16 out of 60 conformations are /3-hairpin structures. The results 
confirm that the original AMBER parm94 strongly favors a-helix structures and that the original AMBER parm96 
slightly favors /3-hairpin structures (because this peptide should form a- helices). PS However, for both AMBER 
parm94 and AMBER parm96 modified to enhance a-helix-forming tendency, almost all (60 and 59) conformations 
are a-helix structures and there are no conformations with /3-hairpin structures. For AMBER parm94 and AMBER 
parm96 modified to enhance /3-sheet-forming tendency, on the other hand, about half the conformations exhibit 
/3-hairpin structures and no a-helix structures are found. 

In Fig. 01 we show five (out of 60) lowest-energy final conformations of C-peptide obtained by the simulated annealing 
MD simulations for the six cases. According to Tabic [V] there are almost 100 % a-helix formations in three cases, 
namely, the original AMBER parm94, the AMBER parm94 modified towards a-helix, and the AMBER parm96 
modified towards a-helix. Here in Fig. 0] we see the differences in the three cases: the modified AMBER parm94 and 
parm96 favor a-helix structures more than the original AMBER parm94 in the sense that the obtained helices are 
more extended (and almost entirely helical) in the former cases. Moreover, we see clear /3-hairpin formations with 
extended /3 strands in the cases of the AMBER parm94 and AMBER parm96 modified towards /3-sheet. 

In Table IVII the numbers of obtained conformations (final conformations) with a-helix structures and /3-hairpin 
structures are listed for the case of G-peptide, which is known to form /3-hairpin structures by experiments. The 
results are similar to those in Table [V] We see that all 60 conformations are a-helix structures in the cases of the 
original AMBER parm94 and the AMBER parm94 and AMBER parm96 that were modified to enhance a-helix 
structures. For AMBER parm94 and AMBER parm96 that were modified to enhance /3-sheet-forming tendency, the 
conformations exhibit /3-hairpin structures and no a-helix structures are found. One difference in the results of the 
two Tables is that the original AMBER parm96 clearly favors /3-hairpin structures in the case of G-peptide, while 
in the case of C-peptide about the same tendency is observed for a-helix and /3-hairpin with this original force field. 
Our modifications of the force fields resulted in the expected changes in the secondary-structure formations. This is 
again more clearly shown in Fig. [S] where five lowest-energy final conformations of G-peptide are displayed for the six 
cases. Overall features are the same as in Fig. 0] 

Therefore, regardless of the secondary-structure-forming tendencies of the original force fields, our modifications of 
the backbone-torsion-energy term succeeded in enhancing the desired secondary structures. 



IV. CONCLUSIONS 



In this article, we proposed to use a double Fourier series in two backbone dihedral angles, <f> and ip, for the 
backbone-torsion-energy term in protein force fields. This is a natural generalization of the conventional torsion- 
energy terms. It is particularly useful in controlling secondary-structure-forming tendencies, because any function 
in the Ramachandran space can be expanded by this double Fourier series. We can easily modify the existing force 
fileds so that specified secondary-struture-forming tendencies are enhanced, by lowering the backbone-torsion-energy 
surface in the corresponding region of the Ramachandran space. We demonstrated this taking the examples of 
AMBER parm94 and AMBER parm96 force fields. 

Besides the above "manual" adjustment of the force fields, we can also apply our force-field refinement method 
[T2I liH , Il4| to this double Fourier series. Namely, we can determine the values of the Fourier coefficients so that the 
forces acting on the atoms with the coodinates of the PDB database become minimal. Work is underway in this 



direction. 
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TABLE I: Torsion-energy parameters for the backbone dihedral angles cj> and ip for AMBER parm94 and AMBER parm96 in 
Eq. @. 







m 


v m 


(kcal/mol) -y m 


(radians) 


n 


V n 
? 


(kcal/mol) 7„ 


(radians) 


parm94 


2 




0.2 


TV 


1 




0.75 


TV 












2 




1.35 


TV 












4 




0.4 


7T 


parm96 


1 




0.85 





1 




0.85 







2 




0.3 


TV 


2 




0.3 


7T 



TABLE II: Fourier coefficients c, normalization constants a, and the basis functions x((f>, ip) for the double Fourier series of the 
backbone torsion energy £(<j>,ip) in Eqs. and JHJ. 



c 


a 




x(cf>, tl>) 


a 


4tt 2 




1 


6m 


27T 2 




cos m<f> 


Crn 


2TT 2 




sin m(j> 


d„ 


27T 2 




cos 


e n 


2TT 2 




sin m/; 


fmn 


IT 2 


cos rri(i 


i cos 


Qjnn 


IT 2 


cos rm 


i sin nijj 


hmn 


7T 2 


sin mq 


i cos nV> 




IT 2 


sinmi 


/> sin nip 



TABLE III: Fourier coefficients in Eq. 1)12^ obtained from the numerical evaluations of the integrals in Eq. 0. "org94" stands 
for the original AMBER parm94 force field. "mod94(«)" and "mod94(/3)" stand for AMBER parm94 force fields that were 
modified to enhance a- helix structures and /3-sheet structures, respectively, by Eqs. I10H and IIH . The bin size e is the length 
of the sides of each unit square cell for the numerical integration in Eq. JHJ . 



bin size e 




1° 










10° 






coefficient 


org94 


mod94(a) 


mod94(/3) 


org94 


mod94(a) 


mod94(/3) 


a 


2.700000 


2.308359 


2 


.308359 


2.700000 


2. 


.308370 


2 


.308371 


hi 


0.000000 


-0.330937 





.390575 


0.000000 


-0. 


.331053 





.390521 


bz 


-0.200000 


-0.101549 


-0 


.157968 


-0.200000 


-0 


.101513 


-0 


.157985 


Cl 


0.000000 


0.509599 





.465469 


0.000000 


0. 


509517 





.465404 


C2 


0.000000 


0.221123 


-0 


.238372 


0.000000 


0. 


.221100 


-0 


.238279 




-0.750000 


-1.164401 


-0 


.401480 


-0.750000 


-1. 


.164500 


-0 


.401437 


d 2 


-1.350000 


-1.333115 


-1 


.267214 


-1.350000 


-1. 


.333073 


-1 


.267170 




0.000000 


0.444390 


-0 


.497739 


0.000000 


0. 


.444289 


-0 


.497800 


e2 


0.000000 


0.241460 





.227452 


0.000000 


0. 


.241451 





.227573 


Ai 


0.000000 


-0.342789 


-0 


.340247 


0.000000 


-0. 


.343087 


-0 


.340249 


gu 


0.000000 


0.367596 





.485922 


0.000000 


0. 


367697 





.485925 




0.000000 


0.527849 


-0 


.405490 


0.000000 


0. 


.527949 


-0 


.405492 


in 


0.000000 


-0.566049 





.579100 


0.000000 


-0. 


.565751 





.579103 
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TABLE IV: Fourier coefficients in Eq. 11211 obtained from the numerical evaluations of the integrals in Eq. ©. "org96" stands 
for the original AMBER parm96 force field. "mod96(a)" and "mod96(/3)" stand for AMBER parm96 force fields that were 
modified to enhance a-helix structures and /3-sheet structures, respectively, by Eqs. H1UH and (I1H . See also the caption of 
Table ITTT1 



bin size e 






1° 












10° 






coefficient 




org96 


mod96(a) 


mod96(/3) 




org96 


mod96(a) 


mod96(/3) 


a 


2 


.300000 


1.908359 


1. 


.908359 


2. 


.300000 


1. 


.908370 


1 


.908371 


hi 





.850000 


0.519063 


1. 


.240575 





.850000 


0. 


.518947 


1 


.240521 


62 


-0 


.300000 


-0.201549 


-0 


.257968 


-0 


.300000 


-0 


.201513 


-0 


.257985 


Cl 





.000000 


0.509599 


0. 


.465469 





.000000 


0. 


509517 





.465404 


C2 





.000000 


0.221123 


-0 


.238372 





.000000 


0. 


.221100 


-0 


.238279 


d x 





.850000 


0.435599 


1. 


198520 





.850000 


0. 


435500 


1 


.198563 


ill 


-0 


.300000 


-0.283115 


-0. 


.217214 


-0 


.300000 


-0. 


.283073 


-0 


.217170 


ex 





.000000 


0.444390 


-0. 


.497739 





.000000 


0. 


.444289 


-0 


.497800 


ei 





.000000 


0.241460 


0. 


.227452 





.000000 





.241451 





.227573 


fxx 





.000000 


-0.342789 


-0. 


.340247 





.000000 


-0. 


.343087 


-0 


.340249 


gxx 





.000000 


0.367596 


0. 


.485922 





.000000 


0. 


.367697 





.485925 


/'ii 





.000000 


0.527849 


-0. 


.405490 





.000000 


0. 


.527949 


-0 


.405492 


in 





.000000 


-0.566049 


0. 


.579100 


0. 


.000000 


-0. 


.565751 





.579103 



TABLE V: Number of final conformations with secondary structures obtained from the folding simulations of C-peptide. The 
total number of folding simulations in each case was 60. See the captions of Tables ITTT1 and ITVI "a-helix" stands for the number 
of conformations including the amino acids which were identified to be "H" (= a-helix) by DSSP. "/3-hairpin" stands for the 
number of conformations including the amino acids which were identified to be "B"(= residue in isolated /3-bridge) or "E"(= 
extended strand, participating in /3-ladder) by DSSP. 



secondary structure org94 mod94(a) mod94(/3) org96 mod96(a) mod96(/3) 
a-helix 60/60 60/60 0/60 14/60 59/60 0/60 
/3-hairpin 0/60 0/60 29/60 16/60 0/60 34/60 



TABLE VI: Number of final conformations with secondary structures obtained from the folding simulations of G-peptide. See 
the caption of Tabic IVl 



secondary structure 


org94 


mod94(a) 


mod94(/3) 


org96 


mod96(a) 


mod96(/3) 


a-helix 


60/60 


60/60 


0/60 


12/60 


60/60 


0/60 


/3-hairpin 


0/60 


0/60 


31/60 


23/60 


0/60 


22/60 
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FIG. 1: Backbone-torsion-energy surfaces of AMBER force fields. The backbone dihedral angles <j> and if) are in degrees, (a) 
and (b) are those of the original AMBER parm94 and the original AMBER parm96, respectively, (a') and (b') are those of (a) 
and (b), respectively, that are expressed by the truncated double Fourier series in Eq. (1121 . The contour lines are drawn every 
0.5 kcal/mol. 
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FIG. 2: Backbone-torsion-energy surfaces of AMBER force fields that were modified by Eqs. II ( )li and pip , (al) and (a2) are 
those of AMBER parm94 force fields that were modified to enhance a-helix structures and /?-sheet structures, respectively, 
(bl) and (b2) are those of AMBER parm96 force fields that were modified to enhance a-helix structures and /3-sheet structures, 
respectively. See also the caption of Fig. Q 
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FIG. 3: Backbone-torsion-energy surfaces of AMBER force fields that were modified by Eqs. 1101 and 11 U and expanded by 
the truncated double Fourier series in Eq. ill '21 . (af) and (a2') are those of AMBER parm94 force fields that were modified to 
enhance a-helix structures and /3-sheet structures, respectively. (bV) and (b2') are those of AMBER parm96 force fields that 
were modified to enhance a-helix structures and /3-sheet structures, respectively. See also the caption of Fig. Q 
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FIG. 4: Five lowest-energy final conformations of C-peptide obtained from six sets of 60 simulated annealing MD runs, (a) 
and (b) are the results of the original AMBER parm94 and the original AMBER parm96 force fields, respectively, (al') and 
(a2') are those of the truncated double Fourier series of AMBER parm94 force fields that were modified to enhance a-helix 
structures and /3-sheet structures, respectively, (bl') and (b2') are those of the truncated double Fourier series of AMBER 
parm96 force fields that were modified to enhance a-helix structures and /3-sheet structures, respectively. The conformations 
are ordered in the increasing order of energy for each case. The figures were created with Molscript |3^1 and Raster3D|34|. 
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FIG. 5: Five lowest-energy final conformations of G-peptide obtained from six sets of 60 simulated annealing MD runs, (a) 
and (b) are the results of the original AMBER parm94 and the original AMBER parm96 force fields, respectively, (al') and 
(a2') are those of the truncated double Fourier series of AMBER parm94 force fields that were modified to enhance a-helix 
structures and /3-sheet structures, respectively, (bl') and (b2') are those of the truncated double Fourier series of AMBER 
parm96 force fields that were modified to enhance a-helix structures and /3-sheet structures, respectively. The conformations 
are ordered in the increasing order of energy for each case. The figures were created with Molscript |3^1 and Raster3D 34]. 



